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Abstract 



This paper offers a new class of models of the term structure of interest rates. We allow each 
instantaneous forward rate to be driven by a different stochastic shock, constrained in such a 
way as to keep the forward rate curve continuous. We term the process followed by the shocks to 
the forward curve "stochastic strings" , and construct them as the solution to stochastic partial 
differential equations, that allow us to offer a variety of interesting parametrizations. The models 
can produce, with parsimony, any sort of correlation pattern among forward rates of different 
maturities. This feature makes the models consistent with any panel dataset of bond prices, not 
requiring the addition of error terms in econometric models. Interest rate options can easily be 
priced by simulation. However, options can only be perfectly hedged by trading in bonds of all 
maturities available. 
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1 Introduction 



In this paper, we develop a new class of bond pricing models that greatly extends the framework 
of Heath, Jarrow and Morton (HJM, 1992). Our model is as parsimonious and tractable as the 
traditional HJM model, but is capable of generating a much richer class of dynamics and shapes 
of the forward rate curve. Our main innovation consists in having each instantaneous forward 
rate driven by its own shock, while constraining these shocks in such a way as to keep the curve 
continuous. This shock to the curve is termed a "stochastic string", following the physical analogy 
of a string whose shape changes stochastically through time. 

Using stochastic strings as noise source for the dynamics of the forward curve has important eco- 
nomic advantages over previous approaches. Existing term structure models have the same set of 
shocks affect all forward rates. This feature constrains the correlations between bond prices, and, 
therefore, the set of admissible shapes and dynamics of the yield curve. In term structure models 
with state variables (e.g. Vasicek, 1977 and Cox, Ingersoll and Ross, CIR, 1985), having complex 
shapes of the forward curve and corresponding elaborate dynamics requires the introduction of a 
large number of state variables. This makes the models very unparsimonious and virtually im- 
possible to estimate. A consequence of the same problem is the difficulty that these models have 
in fitting the term structure at any given point in time, for given parameters estimated from a 
previous time series of data. 

Similar problems affect the traditional HJM model, where fitting closely the initial term structure 
imposes strong constraints on the dynamics of the forward rates, that are in general not verified as 
the model is successively fitted through time. 

In contrast to these models, in our approach, any (finite) set of bond prices is imperfectly correlated. 
The model is thus fully compatible with any given panel dataset of bond prices, or, to put it 
differently, any (finite) set of bond prices observed at some (finite) sampling frequency is consistent 
with the model. There is thus no need to add observation noise when estimating the model. 

Another interesting characteristic of our approach is that, in general, it is necessary to use a 
portfolio with an infinite number of bonds to replicate interest rate contingent claims. However 
pricing remains simple: interest rate options can in general be priced by simulation and, in some 
cases, in closed form. As with HJM, our approach does not allow the formulation of a partial 
differential equation to price derivatives. 

We provide a detailed treatment of stochastic strings, and their stochastic calculus. The framework 
we use is that of stochastic partial differential equations (SPDE's) and our main tool is the calculus 
of Dirac distributions. Methodologically, we attempt to present the results and their derivations in 
the simplest and most intuitive way, rather than emphasize mathematical rigor. 

Kennedy (1997) and, in recent independent work, Goldstein (1997), propose a similar approch 
to modeling forward rates. Kennedy (1997) simply models the forward rate curve as a Gaussian 
random field. Goldstein (1997) uses a model similar to ours, letting the forward rate curve be 
shocked by two of the strings that we analyse as special examples of our framework. 

The paper is organized as follows. In the next section we define the primitives of our approach 
and solve a simple model where forward rates are driven by Brownian motion, to compare with the 
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models driven by stochastic strings. Section 3 defines string shocks and their properties and offers 
the general no- arbitrage condition for forward rate dynamics driven by stochastic strings. Section 
4 discusses te construction of stochastic strings as solutions of SPDE's. In section 5, we present a 
collection of examples of stochastic strings, that present interesting properties. Section 6 presents a 
discussion of the emprical implications of the model, in comparison with traditional factor models. 
In section 7, we show how interest rate derivatives can be priced and hedged in the model. Section 
8 concludes with some directions for future work. 



2 The Traditional Model: Brownian Motion as Noise Source 

We postulate the existence of a stochastic discount factor (SDF) that prices all assets in this 
economy, and denote it by This SDF can be thought of as the nominal, intertemporal, 

marginal rate of substitution for consumption of a representative agent in an exchange economy. It 
is well known that assuming that no dynamic arbitrage trading strategies can be implemented by 
trading in the financial securities issued in the economy is roughly equivalent to the existence of a 
strictly positive SDF. For no arbitrage opportunities to exist, and under an adequate definition of 
the space of admissible trading strategies, the product MV must be a martingale, where V is the 
value process of any admissible self-financing trading strategy implemented by trading on financial 
securities. Then, 



V(t) = E t 



(1) 



M{t). 

where s is a future date and E( [x] denotes the mathematical expectation of x taken at time t. In 
particular, we require that a bank account and zero-coupon discount bonds of all maturities satisfy 
this condition. 

A security is referred to as a (floating-rate) bank account, if it is "locally riskless" .[] Thus, the value 
at time t, of an initial investment of -B(O) units in the bank account that is continuously reinvested, 
is given by 

B(t) = £(0)expj^ r{s)ds\ , (2) 
where r(t) is the instantaneous nominal interest rate. 

We further assume that at any time t riskless discount bonds of all maturity dates s trade in this 
economy and let P(t, s) denote the time t price of the s maturity bond. We require that P(s, s) = 1, 
that P(t,s) > and that dP(t,s)/ds exists. 

Instantaneous forward rates at time t for all times-to-maturity x > 0, f(t,x), are defined by 

/(!■«) -- ak » J %« + ' ) . (3) 



lr This process is also termed the pricing kernel, the pricing operator, or the state price density. We use these terms 
interchangeably. See Dufrie (1996) for the theory behind SDF's. 

2 A security is "locally riskless" if, over an instantaneous time interval, its value varies deterministically. It may 
still be random, but there is no Brownian term in its dynamics. 
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which is the rate that can be contracted at time t for instantaneous borrowing or lending at time 
t + x. We require that the initial forward curve f(0,x), for all x, be continuous. 

Equivalently, from the knowledge of the instantaneous forward rates for all times-to-maturity be- 
tween and time s — t, the price at time t of a bond with maturity s can be obtained by 

P(M) = exp j-^ S * dxf(t,x)X . (4) 

Forward rates thus fully represent the information in the prices of all zero-coupon bonds. 

The spot interest rate at time t, r(t), is the instantaneous forward rate at time t with time-to- 
maturity 0, 

r(t) = f(t,0). (5) 

We use forward rates with fixed time-to-maturity rather than fixed maturity date!^ Modelling 
forward rates with fixed time-to-maturity is more natural for thinking of the dynamics of the entire 
forward curve as the shape of a string evolving in time. In contrast, in HJM, forward rate processes 
disappear as time reaches their maturities. Note, however, that we still impose the martingale 
condition on bonds with fixed maturity date, since these are the financial instruments that are 
actually traded. 

We now introduce a simple model of forward rate dynamics, which is a variation of HJM introduced 
by Musiela (1993) and Brace and Musiela (1994), both for completeness and as a benchmark against 
which to check the models in later sections. 

Let us first simply model all forward rates by 

df(t, x) = a(t, x)dt + a(t, x)dW(t) . (6) 

or, in integral form, 

f(t, x) = /(0, x) + I dv a(v, x) + I dW{v) a(v, x) (7) 
Jo Jo 

From the definition of bond prices, we get at fixed s that 

dlog P(t, s) = f(t, s-t)dt- [' *dy df(t, y) . (8) 
Jo 

3 The model of HJM starts from processes for forward rates with a fixed maturity date. This is different from what 
we do. If we use a "hat" to denote the forward rates modeled by HJM, 

f(t,s) =f(t,s-t) 

or, equivalently, 

f{t,x) = f(t,t + x) 

for fixed s. Musiela (1993) and Brace and Musiela (1994) define forward rates in the same fashion. Miltersen, 
Sandmann and Sondermann (1997), and Brace, Gatarek and Musiela (1995) use definitions of forward rates similar 
to ours, albeit for non-instantaneous forward rates. 
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With Ito's Lemma, we can calculate the dynamics of the prices of bonds with fixed maturity date 



dP(t,s) 
P(t,s) 



f(t,s-t) 



s-t 



dy a(t,y) + - 



o 



s-t ^ 2 

dy a(t, y) 



dt 



Finally, we postulate the following dynamics for the SDFQ 

dM(t) 



M{t) 



-r(t) dt + <j>(t) dW(t) . 



dy a(t,y) ) dW(t) 
(9) 

(10) 



The drift of M is justified from the well-known martingale condition on the product of the bank 
account with the SDF. The process <f> denotes the market price of risk, as measured by the covariance 
of asset returns with the SDF. 

The no-arbitrage condition for buying and holding bonds implies that PM be a martingale in time, 
for any bond price P. Technically, this amounts to imposing that the drift of PM be zero, 

-r(t) + f(t,s-t)- J dy a (t,y) + -(j dy a(t,y)\ - <j){t) (j dy a{t, y)j = . (11) 
Equivalently, for all t and s, with x = s — t 



f(t,x)=r(t)+ / dy a(t,y) 



dya{t,y)) + <j>(t) [ / dy a{t,y) 



(12) 



We can differentiate this no- arbitrage condition with respect to x, and obtain 

a( *' x) = Tx x) + \ {Jo dV V) ) 2 ~ ^ {Jo ^ y) ) } ' (13) 



or, 



a(t, x) = + a(t, x) ( jT dy a(t, y) - 0(t)) . (14) 

where 9 ^q^ is the slope of the forward curve, at time t for time-to-maturity x. This no-arbitrage 
condition clearly shows that the diffusion function a must go to zero when the maturity date goes to 
infinity to ensure the finiteness of aj^| The term 9 ^q^ stems from the parametrization in terms of 
the time-to-maturity. It would be absent in the usual HJM fixed time-of-maturity parametrization. 
Let us stress that the two formulations are completely equivalent. 

We thus get the following arbitrage-free model of forward rate dynamics 

d t f{t, x) = ( d^ g) + a(t, x)\ dt + a(t, x) dW(t) , (15) 

4 M can have other forms than (^). The extension to several Brownian motion terms is trivial, along the lines of 
HJM. More importantly, the SDF can have jumps - that is a Poisson process term - due to discrete changes in the 
arrival of information. 

5 See Dybvig, Ingersoll and Ross (1996) and Jeffrey (1997). 
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denoting 

a{t, x) = a(t, x) (j* dy a(t, y) - <P{t)^j (16) 

This condition can also be found in Musiela (1993) and Brace and Musiela (1994). 

We see that, at some date t, estimating the drift and volatility functions of any forward rate, along 
with the slope of the forward rate curve at that point, recovers the market price of risk process <p^ 
The use of a SDF for pricing is just a useful presentation device. The no- arbitrage condition for a 
single forward rate allows us to construct the dynamics of all forward rates from the knowledge of 
their volatilities. 

In appendix A, we solve equation (|i~5| ) by first rewritting it as 

2M.5^ = *, )+ 4,)fi), (it, 

where dW Jj® is such that, when integrated over time from to t, we obtain Jq dv dW J^ = Jg dW(v). 
We show that the solution of (|l7|) is 

f{t,x) = f{0,t + x) + / dv a(v, t + x - v) + / dW(v) a{v, t + x - v) . (18) 
Jo Jo 

Notice that if a(t,x) = (no fluctuations), we recover the natural no-arbitrage condition f(t,x) = 
f(0,t + x) stating that the instantaneous forward rate is the same at time and at time t for the 
same maturity date t + x. 



3 A Better Model: Stochastic String Shocks 

In this section we develop a more general model of the forward rate curve driven by a stochastic 
string. Instead of the traditional model (^) , we now model the dynamics of the forward rates by 

d t f(t,x) = a(t,x)dt + a(t,x)d t Z(t,x) , (19) 

or, in integral form, by 

f(t, x) = f(0,x) + / dva(v,x)+ / d v Z(v, x) a(v, x) . (20) 
Jo Jo 

where the stochastic process Z(t, x) generalizes to two dimensions the previous one-dimensional 
Brownian motion W(t). The important innovation in ( |i~9| ) is that the stochastic process depends 
not only on time t but also on time-to-maturity x. The notation dtZ(t, x) denotes a stochastic 
perturbation to the forward rate curve at time t, with different magnitudes for forward rates with 
different times-to-maturity. It is straightforward to extend the model to include more than one 
string shock or to combine string shocks with Brownian motion shocks. 

6 This extends naturally to the case where there is a finite number of Brownian increments in the SDF and the 
forward rate processes. In this case, we would need information about the drifts and volatilities of as many forward 
rates as Brownian increments in the model. 

7 In the equation, dW d f^ is an alternative notation for the noise term rj defined later in (^) with covariance (^j|). 
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We stress that (19) is not the infinite-dimensional generalization of the multi- factor HJM model in 
which all forward rates are subjected to the same (possibly infinite-dimensional) set of stochastic 
processes. We have one stochastic process per time-to-maturity. Our process is infinite dimensional 
in the sense that the set of times-to-maturity in the continuous description of the yield curve has 
the order of infinity of the continuum. 

We impose several requirements on Z to qualify as a string shock to the forward rate curve : 



1. Z(t,x) is continuous in x at all times t; 

2. Z(t, x) is continuous in t for all x; 

3. The string is a martingale in time t, E [dfZ(t, x)] = 0, for all x; 

4. The variance of the increments, Var [dtZ(t, x)], does not depend on t or x; 

5. The correlation of the increments, Corr [dtZ(t, x), dfZ(t, x')], does not depend on t. 



We will see that the first two conditions are automatically satisfied by taking Z to be the solution 
of a SPDE with at least one partial derivative in x and t. Condition 3 captures the unforecastable 
aspect of shocks. Condition 4 corresponds to making all shocks affecting the forward curve have 
the same intensity. This intensity can then be changed, or "modulated", by the function a(t,x) 
taking different values at different times to maturity x. The final condition imposes that the 
correlation between shocks two forward rates of different maturities depend only on these maturities. 
Conditions 3, 4 and 5 together make the strings Markovian. We will see that all string shocks 
produced as solutions of SPDE's have Gaussian distributions, which are completely characterized 
by the first two moments. 

To have each forward rate, for each time-to-maturity, be driven by its own noise process, we could 
simply use infinite-dimensional Brownian motion, using a simple Brownian motion as a shock to 
each forward rate, independent from the Brownian motion used to perturb forward rates with 
different times-to-maturity. However, such a model makes the forward rate curve discontinuous as 
a function of time-to- maturity, so that, over time, any two forward rates can become very distant 
from each other. Although a discontinuous forward curve violates no arbitrage conditions, it is 
intuitively unlikely. 

We develop the model, assuming the existence of a string shock that satisfies our requirements. 
In the next sections, we show how such string shocks can be constructed and examine several 
parametric examples. 

We now take the process for the SDF given by 

^) = _ r( £) dt + J" d y y) dtZ (t, y) (21) 

where the pricing kernel is driven by a "weighted sum" of the shocks that affect the forward rate 
curve, and the market prices of risk can in principle be different for each shock to the curve. 
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We follow the steps of the derivation of the no-arbitrage condition on the drift of forward rates of 
section |2|. We have 

rs-t 



y(t,s) = d t logP(t,s) 



f(t,s- t) - a(t,y) dy 
o 



,11 - j 'dy a(t,y) d t Z{t,y) . (22) 
'o 



We need the expression of dP(t, s)/P(t, s) which is obtained from ( P2|) using Ito's calculus. In order 
to get Ito's term in the drift, recall that it results from the fact that, with y stochastic, 



dF ld 2 F 

dtF{y) = — dt y + -—V ar[dty] 



(23) 



so that the differential dfF up to order dt contains contributions from all the covariances of the 
form Cov [dtZ(t, x), dtZ{t, x')] appearing in Var [dty(t, s)]. 



Taking the expectation of the square of the stochastic term in the r.h.s. of (|22|), bond price dynamics 
can be written as 



dP(t,s) 
P(t,s) 



•s—t ^ rs—t rs — t 

f(t,s-t)- I a(t,y)dy + - dy dy'c(t,y,y')a(t,y)a(t,y') 

Z Jo Jo 



dt 



s-t 



where 



dy a(t, y) d t Z(t, y) 



c(t,y,y') = Con [d t Z(t,y),d t Z(t,y')] . 



(24) 



(25) 



and we make use of Var [dtZ(t, x)] = dt. In the next sections, we offer several examples of stochastic 
strings and the correlation functions of their time increments. 

The no-arbitrage condition becomes 

rx ^ rx rx 

-r(t) + f(t,x) - dya(t,y) + - dy dy' c(t,y,y')a(t,y)a(t,y') 
Jo 2 Jo Jo 

(26) 



(27) 



dy / dy' c(t,y,y')(j)(t,y)a(t,y') = 
Jo 

that we can differentiate with respect to x^\ to obtain 

/, v df(t,x 
a(t, x) = 



dx +cr(t,x)ij^ dy c(t,x,y)a(t,y) - dy c(t,x,y)(/)(t,y) 



If c(t, x, y) in (|27| ) is one, we recover the traditional model of a single Brownian motion 
We thus get the following arbitrage-free model of forward rate dynamics 

dtf{t,x) - 



j} r + cr(t, x) M dy c(t, x, y)a(t, y) - dy c(t, x, y)(f>{t, y) 



dt 



In principle, differentiating (g6|) is only warranted if f(t, x) is differentiable in x. This is the case for some strings 
constructed below, such as the integrated O-U sheet, but is not true for others, such as the O-U sheet. In these cases, 
d-^I*'^ is not a function but is similar to an increment of a Brownian motion. However, the results below still hold, 

OX ' ' 

by reading 9 g^ as a convenient notation for the formal integral calculus used. 
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+ a(t, x)d t Z(t, x) 



(28) 



Expression ( |2"8| ) can be solved in a similar way to section |2|. We denote 

A(t,x) = a(t,x) (/ dy c(t,x,y)a(t,y) - dy c(t,x,y)cj>(t,y) 




(29) 



to finally obtain 



f(t,x) = f(0,t + x) + / dv A(v,t + x - v) + [ d v Z(v,t + x- v) a(v,t + x-v) . (30) 




We see that the correlation function of the time increments of the strings is sufficient (together 
with the risk premia) to characterize arbitrage-free dynamics of the forward curve. 

We can still invert the forward curve to obtain the market prices of risk needed for pricing deriva- 
tives. However, we now need to use the entire forward rate curve rather than a smal number of 
forward rates (equal to the number of Brownian motions driving the curve in the traditional model) . 
Note that if we assume that 4>(t,x) = (f>(t), thus constraining the SDF (p|), the "inversion" of the 
forward rate curve to extract the market prices of risk becomes much easier. 

4 Construction of Stochastic String Shocks 

This section shows how to construct stochastic string shocks that can be used to perturb the forward 
rate curve. 

We start by reviewing the construction of Brownian motion from a probabilistic point of view and 
then as the solution to a stochastic ordinary differential equation (SODE). This allows us to obtain 
a natural generalization of Brownian motion to string motion as the solution of stochastic partial 
differential equations (SPDE's). We then formulate the solutions to these SPDE's in terms of the 
Green function of the equation and find constraints on this function to impose our requirements 
for string shocks. 

4.1 Brownian motion as solution of a SODE 

The usual description of stochastic processes in Finance is based on the Brownian motion W defined 



by 




(31) 



normally distributed with 



E [d t W{t)\ 



: 



(32) 



and 



Var [d t W{t)\ 



dt . 



(33) 
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Alternatively, W can be defined as the solution of the following SODE : 

where r/ is white noise, characterized by the covariance 

Cov [r } (t),r l (t')]=6(t-t') , (35) 
and 5 designates the Dirac function.^] 



The SODE ( pif ) describes a particle at position W which is incessantly subjected to random forces 
leading to random variations rj of its velocity. 



The solution of (|34j ) is formally 



W{t) = / dv r]{v) . (36) 



o 



which shows that dt r]{t) is simply a notation for dtW(t), and, as usual, has rigorous mathematical 
meaning only under the integral representation. From (p6|), we see (by the Central Limit Theorem) 
that W will be Gaussian, with mean zero, and we can calculate easily the covariance 

ft i-t' rt ft' 

Cov \W(t), W{t')] = dv dv' Cov[ry(u), r)(v')] = dv dv' 5{v-v) = (tAt 1 ) , (37) 
Jo Jo Jo Jo 

where (t A t') stand for min(i, t'). 

Finally, we define dtW{t) as the limit of [W(t + 5t) — W(t)] when the small but finite increment of 
time 5t, typical of what we have to deal with in real time series, becomes the infinitesimal dt. Using 
(37), we get Var [W(t + 5t) — W(t)] = 5t which recovers (^) in the infinitesimal time increment 
limit. 

The definition of the Brownian motion as the solution of a SODE is very useful to generalize to 
other processes. Maybe the simplest extension is the well-known Ornstein-Uhlenbeck (O-U) process 
U (t) which can be defined as the solution of the following SODE : 

^- = -KU(t)+ V (t) , (38) 

where k is a positive constant. In addition to the random variations 77 of the velocity, the Brownian 
particle is now subjected to a restoring force tending to bring it back to the origin, i.e. mean 
reversion. The formal solution reads 

U(t) = C dv rj{v) e-^^ . (39) 
Jo 

Its covariance is 

Cov [U(t), U(t')} = i-Ye-*!*-*'! - e ~^ t+t 'A , (40) 



9 The Dirac function is a distribution in the sense of Schwartz's theory of distributions, or generalized functions. 
It is such that, for arbitrary f(x), 2 f(x)5(x — xq) — f(xo) if x\ < xq < X2, ^f(xo) if xi — xo or xo = X2 and zero 
otherwise. S(x) is obtained as the limit of functions when their width goes to zero and their height goes to infinity 

in such a way that their integral remains one. An example is the Gaussian limc-^o —7= — e 2^ = 8{x). 
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which goes to j-e K '* * ' at large times -P^l By adding more terms in (j3§|), more complex stochastic 
processes can be easily generated. 

4.2 Stochastic strings as solutions of SPDE's 

We now generalize the stochastic processes driving the uncertainty in the forward rate curve from 
solutions of SODE's to solutions of SPDE's. In the context of differential calculus, this is the most 
natural and general extension that can be performed.]^] The introduction of partial differential 
equations is called for in order to account for the continuity condition and to deal with the two 
variables, time and time-to-maturity, and their interplay. 

Brownian motion depicts, in physical terms, the motion of a particle subjected to random velocity 
variations. The analogous physical system described by SPDE's is a stochastic string. In our 
application of stochastic strings as the noise source in the dynamics of the forward rate curve, the 
length along the string is the time-to-maturity, and the string configuration gives the amplitude 
of the shocks at a given time, for each time-to-maturity. The set of admissible dynamics of the 
configuration of the string as a function of time depends on the specification of the SPDE. 

In the present paper, we restrict our attention to linear SPDE's, in which the highest derivative is, 
in most cases, second order. This second order derivative has a simple physical interpretation : the 
string is subjected to a tension, like a piano chord, that tends to bring it back to zero deformation. 
This tension forces the "coupling" among different times-to-maturity so that the forward rate 
curve remains continuous. In principle, a more general formulation would consider SPDE's with 
terms of arbitrary derivative orders^] However, it is easy to show that the tension term is the 
dominating restoring force, when present, for deformations of the string (forward rate curve) at 
long "wavelengths" along the time-to-maturity axis. Second order SPDE's are thus generic in the 
sense of a systematic expansion] 1 ^ 

The general form of second-order linear SPDE's reads 

, ,d 2 X(t,x) nlt ,d 2 X(t,x) . ,d 2 X(t,x) J . . dX(t,x) dX{t,x)\ 

(41) 

where X is the stochastic string shock we want to characterize. In the present paper, we restrict 
our discussion to linear equations, where F has the form 

, / w N dX(t,x) dX(t,x)\ „ ,dX(t,x) , ,dX(t,x) , . , . , , 

/ it, X, X(t, X), ) = d (*' X ) Qf^ + e (*' X ) g x + X ) + X ) • 

X X (42) 

10 Using 2 (t At') = t + t' - \t - t'\. 

11 Further extensions will include fractional differential equations and integro-differential equations, including jump 
processes. 

12 Higher order derivatives also have an intuitive physical interpretation. For instance, going up to fourth order 
derivatives in the SPDE correspond to the dynamics of a beam, which has bending elastic modulus tending to restore 
the beam back to zero deformation, even in absence of tension. 

13 There are situations where the tension can be made to vanish (for instance in the presence of a rotational 
symmetry) and then the leading term in the SPDE becomes the fourth order "beam" term. 
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s is the "source" term that will be generally taken to be Gaussian white noise r\ characterized by 
the covariance 

Cov [r]{t, x), v(t',x')] = 5(t - t') 8{x - x') , (43) 

where 5 denotes, as before, the Dirac distribution. Expression ( f4l| ) together with (42) is the most 
general linear second-order SPDE in two variables. The solution X exists and its uniqueness is 
warranted once "boundary" conditions are given, such as, for instance, the initial value of the string 
X(0,x), for all x, as well as any constraints on the particular form of equation (f4l|).P| 



4.3 Green function formulation 

All solutions of the above linear SPDE can be characterized by 

r-t r-oo 

X(t,x)=X(0,x)+ dv dy G(t,x\v,y) V (v,y) , (44) 

JO J-oo 

where the Green function G contains all the information in the underlying SPDE and thus of the 
coupling between the different times-to-maturity x. 

Consider a single impulsive source term r](v,y) = 5(v — to)5(y — xo). Then, 

X(t,x) = X(0,x) + G{t,x\t ,x ). (45) 

The Green function G thus describes the deterministic evolution in the future of the process X at all 
times-to-maturity due to an impulsive perturbation that occurs at time to a t the time-to-maturity 
xq. It thus embodies all the coupling between different times-to-maturities. 

We have 

ft roo 



rt roo 

Vax[X(t,x)] = dv dy [G(t, x\v, y)} 2 , (46) 

JO J-oo 



and 



rt+St rod 

Cov [dfX(t, x), dtX(t, x')] = lim / dv dy G(t + 5t,x\v,y) G(t + 5t,x'\v,y) 

St^dtJt J~oc 

rt poo 

+ dv dy[G(t + St,x\v,y)-G(t,x\v,y)][G(t + St,x'\v,y)-G(t,x'\v,y)} (47) 

JO J-oo 

However, this general solution does not satisfy all our conditions for a stochastic string shock to the 
forward curve. We want to impose in addition that X(t, x) is a martingale, and that the variance 
of its increments, Var [d t X(t, x)], does not depend on x. In appendix B, we show that strings that 
satisfy our requirements can be written as 

Z(t, x) = Z(0, x) + f dv f j(X) dy h(v, x, y) v (v, y) , (48) 

JO Jxn 



See for instance Morse and Feshbach (1953). 
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where g, j and h are functions to be constrained in order to have the variance of the increments 
independent of x, and the lower bound can be adjusted in order to get rid of convergence problems 
at infinity. This leads to 

d t Z{t, x) = dt f JX dy h(t, x, y)T]{t, y) , (49) 

Jx 

so that 

Mar [d t Z(t,x)]=dt f dy [h(t,x,y)] 2 . (50) 

The additional condition that Var [dtZ(t, x)\ does not depend on x leads to several solutions. 

A first class of solutions is obtained if we assume that h(t, x, y) is independent of y. The integration 
gives then simply j(x)[h(t, x)] 2 . Equating to a constant (equal to 1 without loss of generality), this 
gives 

h(x) = -±= . (51) 
VJ{ X ) 



The corresponding stochastic string is thus defined by 

1 



d t Z(t, x) = dt —i= [ Ji } dy rj(t, y) . (52) 



and the correlation of the increments is 



c(t, x, x) = Corr [d t Z(t, x),d t Z(t, x')] = y , (53) 

if j(x) < j(x'). The role of x and x' in ( |53[ ) are inversed if j(x) > j(x'). The corresponding SPDE 
reads 

An alternative class of solutions can be obtained by noting that 

where is the inverse function of j defined by iy)) = y, verifies the constraint that 

Var [dtZ(t, x)] be independent of x. This corresponds to 

1 ( d 1 \ 1/2 

d t Z(t,x) = dt—J^ dy ly—\j]-\y)^ v (t,y) . (56) 

A generalization of (|55| ) is 

Mx)] 2 = ~(\J]~ 1 (y)T . (57) 



x a dy 

where a > is arbitrary. In fact, even more generally, consider an arbitrary monotonous function 
I. Then, 
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satisfies the constraint that Var [dtZ(t, x)] does not depend on x. This leads to 



d,Z(t,x) = dt -= I </.;,/—/( [;]- J )f<) '/!'■'/' <W 



The correlation of the increments is 



c(t, x, x') = Corr [d t Z(t, x),d t Z(t, x')] = ^ M , (60) 

if j(x) < j(x'). The role of x and x' in (|(]) are inversed if j(x) > j(x'). This provides a slight 
generalization to (|53| ) since a different function appears in the correlation function and in the 
inequality condition on x and x' . 

The corresponding SPDE is 



4.4 Differentiability in time-to-maturity 

All the processes studied above are continuous both in t and x but not differentiable in either t or 
x. Of special interest is the non-differentiability in x,p| 

We would like now to discuss a class of processes which are non-differentiable in time but dif- 
ferentiable in time-to-maturity, while still keeping an infinite set of shock processes, one for each 
time-to-maturity x.f^\ Economics has little to say with respect to the differentiability of forward 
rate curves. Furthermore, this issue cannot be resolved empirically, since forward rates of very close 
maturities cannot be observed. The use of string shocks that produce differentiable curves is thus 
fundamentally a matter of taste. 

An intuitive strategy to obtain a shock differentiable in time-to-maturity is to integrate any of the 
previously defined strings in the x-variablePl 

Y(t,x) = f X dy Z(t,y) . (62) 
Jo 

By definition, 9Y q^ = Z(t,x). Since Z(t,x) is continuous in x, dY Q^ is also continuous in x and 
Y(t, x) is thus differentiable with respect to x. 

Since Z(t,y) satisfies a SPDE, we see that, by replacing Z(t,y) by 9Y q^ in this equation, that 
Y(t, x) is also the solution of a SPDE. 

The correlation of the increments is 

c(t, x, x) = Corr [d t Y(t, x), dtY(t, x)] = x 
15 We show in appendix C how to determine the differentiability of string shocks. 

16 HJM provides an obvious example of a stochastic process which is non-differentiable in time and differentiable 
in x, but this is a degenerate case as the same stochastic process drives all times-to-maturitics. 
17 Smoother strings can be obtained by higher-order integration. 



'Z(y) + \l{x) 



o V l ( x ) J* V l (y) 



(63) 
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for x < x' , assuming that j(x) is monotonous increasing. The corresponding expressions for other 
cases are straightforwardly derived. 

The variance is 



Var [d t Y(t, x)\ = dt x f dy JM- , (64) 
Jo V l \ x > 

which is independent of x if l(x) = cie C2 ( x ~ 1 / x ' , where c\ and C2 are two arbitrary constants. 

5 Parametric Examples 

In this section we examine several examples of stochastic strings that are solutions to SPDE's, and 
consider their interest as a noise source in a model of the dynamics of the forward rate curve. 

5.1 The Brownian sheet 

We start by discussing the simplest example of a SPDE, which however does not satisfy our criteria. 
The example is however useful as it is the basis for the derivation of other string processes. 

The simplest SPDE in the class given by fl4l| ) and fl42| ) is 

d 2 W(t,x) 



dtdx 



r)(t,x), (65) 



where rj(t, x) is a white noise both in time and x, characterized by the covariance function (^). By 
inspection of (|6^) , the first order time derivative ensures that W will be a Brownian motion in time 
at fixed x, while the introduction of the partial derivative with respect to x ensures the continuity 
of W with respect to x. The solution of ( |65[ ) reads 



W(t,x) = dy dv 7](v,y) . (66) 
Jo Jo 

This process (|66| ) is known as the Brownian sheetQ and has the following covariance, which can be 
readily obtained using (^) with (|^) : 

Cov [W(t, x), W(t', x')] = (t A t') (x A x) , (67) 

where (t At') is the minimum of t and t' . 

For the calculus used above in the implementation of the no-arbitrage condition on forward rate 
dynamics, we need the correlation of the time increments of the string. The covariance is the limit 
when 5t — ► dt of Cov [W(t + St, x) — W(t, x), W(t + St, x') — W(t, x')\. Its calculation reduces to 
that of four terms of the form (|67| ) that simplify to 



Cov [d t W(t, x),d t W(t, x')] = (x A x) dt . (68) 



3 See Walsh (1986). 
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or, stated in terms of the correlation, 

c(t, x, x) = Corr [d t W(t, x),d t W(t, x')] = (x A x') 



(69) 



The property that the variance of W(t, x) is not homogeneous in x, makes this process unsatisfactory 
as a stochastic perturbation in (|l9|): the instantaneous spot rate (x = 0) would have no volatility 
in this model since W(t,0) = 0, and the variance of forward rates would increase linearly with x. 



A simple modification of the Brownian sheet can be obtained by making j{x) = x in (|52|), leading 
to 

d t W(t,x) =dt-= [ dy rj(t,y) . (70) 
Vx Jo 



This process is just like the Brownian sheet, however rescaled by the factor to ensure homo- 
geneity in x. This form could have been guessed directly from the expression ([37]) of the correlation 
of the Brownian sheet. The corresponding SPDE is 

d ( ^dW(t,x)\ . . 

1 G ^-t)= V (t,x) , (71) 



dx V dt 

and 



Corr 

for x < x' and for x > x' . 



d t W{t,x),ckW(t,x') = J- , (72) 



5.2 The Ornstein-Uhlenbeck sheet 

A natural model for Z(t, x) is one similar to the Brownian sheet but in which the variance is 
the same for all times-to-maturity. We look for the string with the same variance for all x and 
covariance between x and x' that depends only on the difference (x — x'). 



The O-U string satisfies ( |52[ ) with j(x) = e 2KX , which leads to 

U(t, x) = e- KX W(t, e 2KX ) = e~ KX f dy C dv rj(v, y) . (73) 

Jo Jo 

The corresponding SPDE is 



which we can simplify to 

d 2 U(t,x) dU(t,x) r— . . 

' k y ' =Vteri(t,x) . (75) 



dxdt dt 
By inspection of (^), we get the following properties 

Cov [U(t, x), U(t', x')] = (t A t') e - K \ x ~ x '\ , (76) 
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and 

c(t, x, x') = Corr [d t U{t, x),d t U(t, x')] = e'^'^ . (77) 

Expressions ( |76| ) and (77) both show that, for fixed x (x = a/), U(t,x) is a Brownian motion in 
time with variance independent of x. This is the property that we were looking for, so that all 
times-to-maturity have a priori the same volatility for the driving stochastic process. □ We also 
learn that the time increments dtU(t,x) and dtU(t,x') are correlated accross times-to-maturity: 
this is necessary to obtain the continuity condition along x. Again, if the variations along x were 
delta-correlated, the process would consist of independent one-dimensional Brownian motions, one 
for each x, thus leading to an almost everywhere discontinuous process. The correlations along x 
are the strongest between close times-to-maturity and vanish exponentially between distant times- 
to-maturity. The parameter k allows us to span a set of models, from simple Brownian motion 
to a stochastic string. Indeed, for k — > 0, all times-to-maturity are so strongly coupled that they 
all consist of the same single dimensional Brownian motion. This limit K — > thus recovers the 
single factor model driving simultaneously all times-to-maturities. Figure 1 shows an example of 
the increment in the O-U sheet for different values of k. We see that, as k increases, the shape of 
the string has less and less humps. 



There is another way to construct the O-U sheet. Let us come back to the definition (|7q) of U(t, x) 
and make the change of variable y — ► e 2KZ in the integral. We get 



U{t,x) = 2 K e~ KX [ X dze 2KZ f dv r,{v,e 2KZ ) . (78) 

J-oo JO 

We can write rj(v, e 2KZ ) = ^-7=fj(t, z), where fj(t, z) is still white noise Cov [rj(t, z),r)(t' , z')] = 5(t — 



2k 

t') 8{z — z'). Dropping the hat on the noise, we get 

rx rt 



U(t, x) = V2^e~ KX r dz e KZ f dv rj(v, z) . (79) 

J-oo JO 

We then write J* = /° + Jq. The first integral gives a contribution 

e~ KX f dv [ dz e KZ V2^ rj(v,z) = e~ KX C dv r}(v) , (80) 

Jo J-oo Jo 

where r)(v) = \/2k J ^ dz e KZ ri(v, z) is a simple Brownian process of unit variance. We have kept 
the same notation rj, but it is clear that it now depends only on time. This is thus the noise process 
that is common to all times-to-maturity. If considered alone, this would recover the single factor 
model driving simultaneously all times-to-maturity. Putting everything together, we obtain the 
alternative formulation of the O-U sheet 

U(t, x) = e~ KX f* dv rj{v) + v 7 ^ e~ KX f dv [* dz e KZ r){v, z) . (81) 
Jo Jo Jo 

This form exemplifies the relationship with the standard one-dimensional O-U process, in the 
sense that it shows the exponential dependence of the process on innovations at different times-to- 
maturity. This is thus the generalization in the time-to-maturity dimension of the classical O-U 

19 Different volatilities for forward rates with different times-to-maturity can be obtained by choice of the weighting 
factor a(t,x) in jl9|). 
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process usually denned in the time direction. The analogy with the one-dimensional O-U process 
is also clear from the SPDE formulation ( |75| ) where it is seen that the exponential dependence is 
reflected in the mean reversion term which plays a role similar to the mean reversion term 
-kU in (H). 

This process, written in the form (|l|), has been used by Kennedy (1997) and Goldstein (1997) 
to model forward rates.f^] Our approach clarifies the special role played by this process and its 
relationship with SPDE's. 



5.3 The integrated Ornstein-Uhlenbeck sheet 

As an interesting example, introduce the "integrated O-U sheet", defined by 

V 1 (t,x)= f dyU(t,y) , (82) 



o 



where U(t,y) is the O-U sheet defined in (|7^). By definition, 9Vi q^ = U(t,x). Since U(t,x) is 
continuous in x, is also continuous in x and Vi(t,x) is thus differentiable with respect to x. 

since U(t,y) satisfies the SPDE ©, we see, by replacing U(t,y) by ^ in this equation, that 

Vi(t, x) is the solution of the following SPDE 

d 3 Ui(t,x) d 2 Vi{t,x) ,— 

This is a SPDE of order three. 
We obtain 

c(t,x,x , ) = CoTT[d t V 1 (t,x),d t V 1 (t,x')] = f dy [* dy' Corr [U(t, y), U(t,y')] 

Jo Jo 

= - (2(x A x') - -(1 - e- KX )(l - e- KX ')) , (84) 

K \ K J 

where we have used (f77|). 

Inspired by the alternative formulation (|8l] ) of the O-U sheet U(t, y) in terms of the exponential 
kernel, we define another integrated O-U sheet V2(t,x) by 

V 2 (t, x) = nV2e~ KX C dy e Ky U(t, y) , (85) 
Jo 

Again by construction, 9V2 q^ x ^ is continuous and thus V2(t, x) is differentiable in x. By differentia- 
tion of (|85|) with respect to x, we get 

= - K v(t,x) + kV2 U{t,x) . (86) 

ox 



20 Kennedy (1997) only considers the case of a constan volatility function for the forward rates, whereas Goldstein 
(1997) allows for a general volatility function. In contrast to our approach, both authors model forward rates of fixed 
maturity date. 
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The SPDE obeyed by V2(t,x) is simply obtained by putting U(t,x) = ~^^r^^r^ + K V2{t,x)} in 
(|75|). This is again a parabolic SPDE of order three. 



The correlation of the time increments of this process looks nicer than (H) : 

c(t,x,x') = Corr [ <i t U 2 (*,£), d*U 2 (t, x')] 

= 2k 2 e- K(x+x,) f dy f X dy' e Kiy+y,) Corr [dtU(t, y), d t U(t,y')] = (1 + k\x - z'De""!*-*'! . (87) 



Figure 2 shows an example of the increment in the integrated O-U sheet for different values of n. 
We see that the string is always smooth and, as k increases, its shape has less humps. 

This process has been independently used by Goldstein (1997) to model forward rates. Our ap- 
proach shows that it is one particular example among an arbitrarily large class of processes exhibit- 
ing similar properties. 



5.4 The term structure of correlations 



Empirical observation shows that correlation between two forward rates, with maturities separated 
by a given interval, increases with maturity. What are the string processes that exhibit this prop- 
erty? Let us start from the parametrization obtained by assuming that h(t, x, y) is independent of 
y, which leads to (|53|). Consider two time-to-maturities x and x' = x + Sx, where 5x > is small 
compared to x. We can then expand j(x') = j(x) + 5x d ^^} + higher order terms. Expanding the 
square root in (jag), we get 

c(t, x, x') = Corr [d t Z(t, x),d t Z(t, x')] « (l - y ^|^) • (88) 

We see that Corr [dtZ(t, x), dfZ(t, x')\ increases with x if decreases (while remaining posi- 

tive) as x increases. An interesting situation is when — > + as x gets large. In this case, 
the increments become perfectly correlated for large time-to-maturities. 

Parametrizing j(x) = e % ^ x \ the above condition implies that i(x) must be a function which increase 
slower than x, in such a way that its derivative decreases with x, that is, i(x) must be concave. In 
order words, j(x) must grow slower than an exponential. This is not the case of the O-U sheet or 
the integrated O-U sheet. However, the previous example (y2\), where j{x) = x, qualifies. In fact, 
any expression like j(x) = e x with an exponent < a < 1 qualifies. This corresponds to 

c(t, x, y) = Corr [d t Z(t, x), d t Z(t, y)] = e — VI . ( 89 ) 

For instance, take j(x) = e 2K ^ x . This corresponds to a correlation function 

c(t, x, y) = Corr [d t Z(t, x),d t Z(t, y)] = e — Iv^-v^l ( 90 ) 

which again works fine. Figure 3 shows the correlations between forward rates at maturities spaced 
by different time intervals, as the maturies increase. We see that correlation increases both with 
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the proximity of the maturities of the forward rates and with the times-to-maturity of the forward 
rates. 

These models make it very easy to fit any covariance matrix of the increments of instantaneous 
forward rates. The function a can be used to fit the variances, and the parameters k in (|90| ) or k 
and a in ( |89[ ) can be used to fit the correlations. 

6 Empirical Implications 

Traditional factor models of the term structure, are only compatible with samples that include 
at most as many bonds as there are factors in the model. In order to be able to use data on 
more bond prices to estimate these models, it is necessary to add error terms to their econometric 
specification. For econometric tractability, these errors must be assumed to be independent from 
the factors p 1 ] This approach is correct if the incompatibility of the data with the model is due 
to observation error, possibly due to bid-ask spreads or nonsynchronous observations. However, if 
the incompatibility is due to a misspecified model, the error terms will not be independent of the 
factors and the econometric model will be inapropriate. 

Models with stochastic string shocks can be estimated without any such error terms. The models 
are compatible with any sample of forward rates (or bond prices) of finitely spaced maturities, 
taken at some finite sampling interval. This is so, because, for any parametric specification, there 
is always a possible path for the string shock over a finite interval that can lead from the forward 
curve at the begining of the interval to the forward curve at the end of the interval. This realization 
may be highly unlikely, but it is always possible. The estimation exercise is thus one of finding the 
most likely parameters, given a set of movements of the forward curve over time. 

7 Option Pricing and Replication 

In our model, derivatives can in general be priced by simulation. We just need to simulate the 
dynamics of the forward curve under the risk adjusted probability measure. Define a new string 
Z* with dynamics 

POO 

d t Z*(t, y) = d t Z(t, y)- dy c(t, x, y)4>(t, y) . (91) 

Jo 

Under the risk adjusted measure, this string is a martingale. Asset prices discounted at the bank ac- 
count are also martingales under this probability measure. Thus, the price at time t of a (European, 
path-independent) contingent claim with payoff <3?(s) at time s > t is 

(92) 

It is sometimes convenient to use the s-maturity bond price as numeraire. In order to change drifts 
21 See Pearson and Sun (1994) or Chen and Scott (1993). 



*(*) 



$(s)e~ 



■ J du r(u) 



19 



for the s-forward risk neutral probability measure, note that 

rs—t rs~t 

v(t,s)= / dx I dy c(t,x,y)a(t,x)a(t,y) (93) 
J o Jo 

is the instantaneous variance of the s-maturity bond price. Then, defining a new string Z s with 
dynamics 

d t Z s (t, y) = d t Z*(t, y) + v(t, s) dt , (94) 
we can price the contingent claim as 

*(t) = P(M)E£[$( a )] . (95) 

Now the expectation is taken with respect to the probability measure under which Z s is a martin- 
gale. 

Unfortunately, we cannot use the Feynman-Kac Theorem to obtain a PDE for pricing contingent 
claims. In our model, there is no finite set of state variables that determine the value of interest 
rate derivatives. Derivatives will in general depend on the full history of (string) shocks to the 
forward rate curve. This shows that to hedge them, it is necessary to use the full set of bonds of 
all maturities. 

Closed form solutions for option prices can be obtained for Gaussian cases, where the volatilities of 
forward rates are deterministic. Kennedy (1994) gives a formula for pricing European call options 
on bonds when the string used is the O-U sheet and a(t, x) = ae~ ax , with a and a positive constants. 



8 Future Work 

In future work, we intend to estimate and test different parametrizations of our model. Another 
avenue of research that we are currently exploring is the modelization of forward rates directly as a 
stochastic string, i.e., as the solution of a SPDE, rather than the approach followed in the present 
paper of using stochastic strings multiplied by volatility functions to shock forward rates. This 
new approach is particularly promising for a parsimonious description of forward rate curves with 
non-linear dynamics. 
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A Solution of Traditional Model 



In this appendix, we obtain the arbitrage-free process for forward rates of the traditional model of 
section ||[ 

We perform a change of variable (t,x) — ► (r = x + t,£ = x — t) and denote /(r, £) = f(t,x). We 
then have 

a/(t,x) _ dT df(r,o ag a/Q-.fl _ a/(r,p a/(r,o 
at at dr at ae #r ae ' 1 j 

and, similarly, 

a/(t,x) = dr dfj^Q | ag a/(r,g) = a/(r,p | a/(r,g) 
ax ax ar dx a£ ar a£ 

Replacing these in the l.h.s. of (pTj), we get 

g/fog) g/M O a/(r,o 

~~flt ax - " ~ 2 ^T~ • (98) 

Equation ( |T7| ) thus becomes 



df(r,0 1 /r-£ r + a 1 /r-e r + n diy(t) 

rCM — = — , — = — — tt(T 



a^ 2V2'2y2V2'2y eft 



(99) 



T — c 



where we have expressed the values of t and x in terms of the new variables r and £ : t - , 
and x = — 1£. This change of variable allows us to get an ODE from the initial PDE and its 
straightforward integration gives 

Now, by definition of the change of variables, /(t, £) = /(t, x). We have a free choice for £o- Since 
we want to express /(t,x) as a function of /(0, .), £o is chosen such that it corresponds to to = 0. 
To see what this leads to, we use the definition /(r, £o) = /(*o = I -^ B ~7 For to to be zero, 
this implies that £o = T an d thus T ~^° = r = x + 1. The l.h.s. of ( |10C| ) thus reads 



/(r, - /(r, Co) = f(t, x) - /(0, x + t) (101) 



To tackle the r.h.s. of ( 100[ ), we perform the change of variable from w to v = Then 



£±2£ = T — t> and J| o J ! J_ =x ._ ) _ i = —2 / dv. Finally, we get 

f(t,x) = /(0, x + t) + / dv a(v,t + x-v)+ [ dW(v) a{v,t + x - v) . (102) 
Jo Jo 

Let us now comment on the validity of this derivation when f(t,x) is not differentiable in x (and 
in t as it is usually). As a teaching example, consider the SPDE 

dW dW 

ST - 17 = • < 103 > 
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Performing the same change of variable as above, we get 

^--itoO. (104, 

where 

r)(T,0 = v((T-0/2,(r + 0m • (105) 

The solution of ( |104| ) reads 

W(j,0 = — dyr)(T,y) = -- dyr ] ((x + t-y)/2,(x + t + y)/2). (106) 
2 Jo 2 Jo 



One readily checks that fllOq) obeys ( |103|) . The terms formally involving the derivatives of the 
noise r\ would take sense when discretizing the equations. 



B Requirements for String Shocks 
B.l Martingale condition for strings 



We want Et[dfX(t, x)\ to be equal to zero. From (f44[), we have 

d t X(t,x) = [*dv dy dG ^ v ^ rj {V: y) + |°° dy G(t,x\t,y) V (t,y) . (107) 



o 



The expectation E t [dtX(t, x)] conditional on the realization of X(t,x) = g(x), where g(x) is a 
specified function, can be written as 

E t [d t X(t,x)] = J... J V{ V }V({ V }) d t X(t,x) s(x{t,x)-g(x)j , (108) 

where J ... J T> {r/} denotes the functional integraQ over all possible realizations of the noise 77, 
V{{rf\) is the probability density function weighting each realization of the noise, and the term 

5^X(t, x) — g(x)^j ensures that the functional integral is carried over all possible ?7's restricted to 

satisfy the specific realization of X(t, x). With delta correlated shocks rj, V({rj}) is Gaussian, so 
that 

p({ v }) = P e ~f dv f dy c ~ 1(v ' y) Mv ' y)]2 , (109) 

where 2C~ 1 (v,y) is the inverse of the variance of rj(v,y) and we allow it to vary with time and 
time-to-maturity. 

Taking the Fourier transform of Et[dtX(t, x)] with respect to g{x) gives 

E t [d t X(t,x,k)} = J -J V{n}V{{n}) dtX(t,x) e ik x ^ , (110) 



22 To get an intuitive representation, imagine that time is discrete. Then, the integrals are carried out over the 
variables r\i for time i = 1 to t. 
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where k is the conjugate of g(x) in the Fourier transform. We notice that the second term in the 



r.h.s. of (107) does not contribute to Et[dtX(t,x,k)] since it contains the innovations rj(t,y) that 
are posterior to those contributing to X(t,x). The gaussian integrals over the centered innovations 
rj(t, y) thus vanish. 



The integrand in (|11C| ) contains the exponential term 

«p(-/*/*[ t r'(,,,) W .,,)]»-*c(M.,y)^,»)" 

C~ l {v,y) [r](v,y)] 2 — ik G(t,x\v,y) rj(v,y) can be factorized as 

C' 1 ^, y) [rj(v, y)] 2 - ik G(t, x\v, y) n(v, y) = 
n -x ( s(( s -u G(t,x\v,y) \ 2 k 2 [G(t,x\v,y)] 2 

C (v,y) ^.y)-^ — — j + , c ^ y) • (HI) 
The Gaussian integrals can be carried over f)(v,y) = i](v,y) — ik fc-i^y) ' wn i cn leads to 

E t [d t X(t,x,k)] = dv dy [ k . - g^l|MjJ e Jfe=i^ . (112) 

Jo J~oc C l {v,y) at 



We look for the condition that the Green function G(t, x\v, y) must satisfy for this expression (112) 
to be identically zero for any t and x,p| 

The condition is simply that 

OGi^yl = Q (U3) 

for all t and x, i.e. G(t, x\v, y) must be independent of t. We note that the dependence in t can still 
appear through dependence in time of possible upper and lower bounds in the integral over y. We 
verify that this condition is observed for the parametric examples of stochastic string shocks given 
in the text. 



The condition ( |113j ) implies that the Fourier transform ( |110| ) of Et[dtX(t,x)] is zero, and thus 
Et[dtX(t,x)] is zero itself. This is a Martingale condition. 

B.2 Markov condition for strings 

An additional requirement (more stringent than the previous one) is that X(t, x) be Markovian, i.e. 
that its increment d±X(t,x) be completely uncorrelated from past time innovations. The general 
form that the stochastic string can take in order to fulfill this condition is 

X(t,x) =X(0,x) + f dv g{v) f ji) 1 dy h(v,x,y) V (v,y) , (114) 

JO J-oo 



If the Fourier transform is zero, the function is zero. 
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where g(v), j(x) and h(v,x,y) are arbitrary functions. This expression obeys the condition ( 113 ) 
and in addition time does not appear in the bounds of the integral over y. As a consequence, from 
dnj), we g et 

d t X(t, x) = dt g(t) [ 3 X dy h{t, x, y) V (t, y) . (115) 



The fact that r](t,y) is immediately posterior to t ensures that the expectation of dtX(t,x) con- 
ditional to the realization of X(t,x) is identically zero. This recovers a result of Kennedy (1997), 
obtained in a different setup. 



B.3 Constraints on the covariance of the increments 

The variance Var [dtX(t, x)\ of the increments must not depend on x. This condition requires 
some additional condition on the functions j(x) and h(v,x,y) that we now derive. Notice that 
Cov [dtX(t,x)dtX(t,x')] does not depend on t if g(t) is constant. 

From (|115p , we get 

Vax[dtX(t,x)]=dt !° ) dy [h(t,x,y)] 2 , (116) 
Jo 

putting g(t) = 1. We have modified the lower bound of the integral in y so as not to have to worry 
about additional constraints in order to ensure convergence at — oo. The condition to impose is 
that ( |116| ) be independent of x, which leads to the solutions given in section 4.3. 



C Differentiability in Time-to-Maturity 

The string processes studied in section 4.3 are continuous both in t and x but not differentiable 
either in t or x. Of special interest is the non-differentiability in x. Consider first the Brownian sheet 
process (|66|) and let us consider 9W q^ = Jq1](v,x). This integral is a continuous sum of random 
variables. The sum of Gaussian random variables is itself Gaussian with variance equal to t. Thus, 
9W Qx' x ^ is not a function but should rather be interpreted as a noise process, discontinuous almost 
everywhere with respect to x. The derivative with respect to x does not exist strictly speaking, but 
only in the sense that its integral is well defined. This non-differentiability holds true for a large 
class of noise processes, not necessarily Gaussian. 

A similar reasoning holds for all the string processes defined in (|4q). 

All the stochastic partial differential equations that we discuss in this paper are linear. For such 
linear equations, a particularly convenient mathematical tool is the Fourier or Laplace transform. 
Let Z(t,x) be one of the stochastic strings. Its Fourier transform with respect to x is 

Z(t, k) = f + °° dx e ikx Z(t, x) . (117) 



The inverse Fourier transform retrieving Z(t, x) from Z(t, k) is 

Z(t, x) = [ +C ° dk e~ ikx Z(t, k) . (118) 
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We have the following useful property : 



dx e 



ikx 



d p Z(t,x) 
dx p 



{-ikf Z{t,k) , 



for an arbitrary positive integer p If r](t, x) is a Gaussian white noise, we have 

Cov [fj(t, k),r)(t', k')] = 5(t - t') 5(k - k') . 



(119) 



(120) 



Consider first the Brownian sheet (65). Taking the Fourier transform of (|65|) leads to 

dW(t,k) i 



-V(t,k) 



whose solution is simply 
We thus get 



Of, 



W(t,k) = W(0,k) + - dTfj(r,k) 
k Jo 



Cov 



5{k - k') 



k 2 



(t A t') 



W(t,k)W(t',k' 

This expression allows us to calculate any desired quantity. For instance, we find 

~dW{t,x) dW(t',x') 



which shows that E 



Cov 



dW(t,x) 
dx 



(t A t') 5(x - x') , 



(121) 

(122) 
(123) 

(124) 



Consider now the O-U process 
process, leads to 



dx dx' 

is infinite, the signature that W(t, x) is not differentiable in x. 

Taking the Fourier transform of the SPDE defining the O-U 
dZ(t,k) 



-ik + k) 



at 



'2k fj(t,k) 



Its solution is 



Z(t,k) = Z(0,k) + 



'2k 



-ik + k Jo 



Following the same steps as before, we get 
'dZ(t,x),dZ(t',x' 



Cov 



dx 



dx' 



(t A 



dr fj(rj, k) . 



dk e lk(x - x,) 



2k k 2 
k 2 + k 2 



(125) 
(126) 

(127) 



The integral still diverges for \k\ — > oo, again expressing the non-differentiability of Z(t,x). 



In contrast, the integrated process (82) which obeys the SPDE (83) has its Fourier transform Y(t, k) 
satisfies 



(k 2 + ikK) 



dY{t,k) 
dt 



r2n fj(t,k) 



(128) 



24 Derivatives are taken in the sense of distributions and we use the usual definition of the Fourier transform of 
distributions. 
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The solution is 



leading to 



Y(t,k) = Y(0,k) 



2k 



dr 77 (r, k) 



Cov 



Y(t,k),Y(t',k') =(tAt 



k 2 + ikn Jo 

2n5(k-k') 



We thus see that 

'dY(t,x) dY(t',x') 



Cov 



2k (t A t') 



k A + k 2 k 2 ' 



(129) 
(130) 



Z^2 1 /.- 2 

— OO ^ \ rv 



vr (t A t 7 ) e~ Klx - x ' 1 , (131) 



<9x cte' 

which is finite, signaling as expected the differentiability of Y(t, x) with respect to x. 
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FIGURE CAPTIONS 



Figure 1 : Sample increment to the O-U sheet, with correlation function ([77]) for different values 
of the parameter k. The time increment is taken to be 1, measured in the same units of time-to- 
maturity. 



Figure 2 : Sample increments to the integrated 0-U sheet, with correlation function (^7j) for 
different values of the parameter k. The time increment is taken to be 1, measured in the same 
units of time-to-maturity. 



Figure 3 : Correlations between forward rates with times-to-maturity separated by different inter- 
vals, for the string shock with correlation function (p0|). The parameter k is taken to be 1. 
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